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We re-analyze the data of the elastic electron proton scattering to look for model- 
independent evidence of the two-photon-exchange (TPE) effect. In contrast to previous 
analyses, TPE effect is parametrized in forms which are free of kinematical-singularity, in 
addition to being consistent with the constraint derived from crossing symmetry and the 
charge conjugation. Moreover, we fix the value of R = Ge/Gm as determined from the 
data of the polarization transfer experiment. We find that, at high Q 2 > 2GeV 2 values, 
the contribution of the TPE effect to the slope of gr vs. e is large and comparable with 
that arising from Ge- It also behaves quasi- linearly in the region of current data, namely, 
in the range of 0.2 < s < 0.95. Hence the fact that the current elastic ep cross section data 
shows little nonlinearity with respect to e can not be used to exclude the presence of the 
TPE effect. More precise data at extreme angles will be crucial for a model-independent 
extraction of the TPE effect. 
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The electric and magnetic form factors Ge and Gm of the nucleon contain essential piece 
of our knowledge of the nucleon structure. Traditionally, Rosenbluth method is used to 
extract Ge and Gm of the proton from the differential cross section of the electron-proton 
{ep) scattering. Within the one-photon-exchange (OPE) approximation, one has 

- g(1 + T) = G\m 2 ) + -Gl{Q% (1) 

Ml Lab TGMott T 

where r = Q 2 /AM 2 and 1/e = 1 + 2(1 + r) tan 2 9^/2, with 9 Lab the laboratory scattering 
angle such that < e < 1. (JMott is the Mott cross section. At a fixed value of Q 2 , the slope 
and the intercept of <tr(£) vs. e is directly related to Ge and Gm, respectively. Through this 
method it is determined that 

ti p R = v P G E (Q 2 )/G M (Q 2 ) « 1, (2) 

where \x v = 2.79 is the proton anomalous magnetic moment. 

Another way to extract the ratio R = Ge/Gm is to measure the polarizations of the recoil 
proton in the elastic ep scattering. Within the OPE approximation, the values of R can be 
extracted according to the following formula [TJ: 




R = -\hhrr<-'!?, (3) 



where Pi and Pt are the polarizations, parallel and perpendicular to the proton momentum 
in the scattering plane, respectively. This method has several advantages over the tradi- 
tional Rosenbluth separation method. The reason is that the ratio of the two simultaneously 
measured polarization components greatly reduces systematic uncertainties. For instance, a 
detailed knowledge of the spectrometer acceptance which plagues the cross section measure- 
ments is not required here. Furthermore, the detailed knowledge of the beam polarization 
and the the polarimeter analyzing power are also not needed because both quantities cancel 
in measuring the ratio of the form factors. Therefore this method has been considered to be 
more reliable. 

Recent polarization transfer experiments at Jefferson Lab [21 [3] give, however, 

fx p R = im p Ge{Q 2 )/G m {Q 2 ) = 1 - 0.13(Q 2 [GeV 2 ] - 0.04), (4) 

which differs substantially from the Rosenbluth results of Eq. ([2]) at high Q 2 . It was suspected 
that there is a systematic uncertainty common to all Rosenbluth data. However, recent 
measurements of high Q 2 ep elastic cross section at the Jefferson Laboratory (JLab) [U [5] 
confirmed the earlier Rosenbluth measurements from SLAC (6j [7] in the same Q 2 region. 

It has been proposed [U [TO] that this mismatch can be explained by the two-photon- 
exchange (TPE) effect not fully accounted for in the standard radiative corrections procedure 
of Mo and Tsai [11]. So far the only observation of possible TPE effect is in the asymmetry 
of scattering of transversely polarized electrons from protons. This observable should vanish 
within the OPE approximation but both SAMPLE experiment at Bates |12] and PVA4 
Collaboration at Mainz [TO] have observed small but non-zero results. This observable is 
related to the imaginary part of two-photon exchange amplitude. On the other hand, the 
unpolarized cross section depends only on the real part of the two-photon exchange amplitude 
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which does not, in general, vanish. It is hence reasonable to expect that the TPE effect would 
contribute to <tr. The crucial question is whether this effect is large enough to explain the 
discrepancy between the results of the Rosenbluth and the polarization transfer methods for 
the ratio R. 

To be more specific, let us take a closer look on how the TPE effect would modify the 
Rosenbluth formula. Including the TPE effect, Eq. ([T]) is rewritten in the following general 
form [8]: 

a R (Q 2 ,e)=G 2 M (Q 2 ) + -G 2 E (Q 2 ) + F(Q 2 ,e), (5) 

r 

where F(Q 2 ,e) is a real function describing the effect of the 17(8)27 interference. Although the 
magnitude of F is suppressed by a em = 1/137 compared to the OPE effect, the e dependence 
of F(Q 2 ,e) is not necessarily small in comparison with eG^/r, particularly at higher Q 2 
region where Ge and 1/r are both small. Hence the slope of the <jr{e) is no longer just 
G 2 e /t but also includes the linear part of the TPE effect. As a result, the validity of Eq. §\§ 
becomes questionable and the extraction of Ge and Gm from cjr(e) needs more sophisticated 
consideration and the values of R = Ge/Gm obtained via the Rosenbluth method should 
also be modified. 

There are currently two ways to study the TPE effect in ur. The first one is to estimate 
F(Q 2 ,e) by calculating the two-photon-exchange ep scattering amplitude with some simple 
hadronic or partonic models (9j [TU1 I14j . However, so far there is no complete calculation 
which is valid at all kinematics. The second way is to analyze the experimental data to 
search for the nonlinearity in ctr with respect to e which is regarded as the evidence of the 
TPE effect I16j since F(Q 2 ,e) is, in general, not a linear function of e. 

All current model calculations [HJ \10\ [T3] show that the existing discrepancy between the 
Rosenbluth and the polarization transfer methods can at least be partially explained by the 
TPE effect [U[T0]. However, several groups [151 116j have claimed that they cannot find any 
conclusive evidence for the the nonlinearity of <jr{e) in their data analyses. The analysis of 
Ref. [TB] set tight limits on the size of the nonlinearity of (Tr(e). Nevertheless, these limits are 
still consistent with the predictions of the model calculations. In fact, this issue was already 
discussed in [8]. It was speculated that there may be a simple explanation of the fact the 
Rosenbluth plot looks linear even if it is strongly affected by the TPE effect. Namely, even 
if the TPE effect in the cross section is not a linear function of e, it is still possible for F to 
behave quasi-linearly within the present experimental uncertainty in the current data region 
(0.2 < e < 0.95). If this is the case, both of the success of Rosenbluth separation in the 
past and the newly discovered discrepancy can be explained because the ur{£) in the certain 
range of e could look like a straight line with its slope a sum of the OPE and TPE effects. 

To explore this possibility, one needs to look for some appropriate parametrization of the 
TPE effect to account for the data of the cross section of the elastic ep scattering. Such 
a parametrization has to possess general characteristics of the TPE mechanism, as derived 
from charge conjugation and crossing symmetry in Ref. [17j . Namely, the last term F in 
Eq. ^ which corresponds to the I7 (g> 27 interference effects should, at a fixed value of Q 2 , 
satisfy the following constraint 

F(Q 2 ,y) = -F(Q 2 ,-y), (6) 
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where y = y - — . In [15], F(Q 2 ,e) is assumed to take the form 

F(Q 2 ,e)~-or -. (7) 

y y 

However, it suffers from the problem that 1/y diverges when y approaches zero, i.e., 9iab - > 0, 
and is hence not acceptable. 

In this letter, we study the TPE effect in an(e) directly through data analysis. We choose 
the data from [7] in order to avoid complication due to different normalizations of different 
experiments. According to the general argument of [8], R = Ge/Gm determined in the 
polarization transfer method is little affected by the TPE effect, at least in the range of Q 2 
which has been determined till now. As a result, it is important to fix the values of R as 
obtained from polarization transfer experiment in any analysis. With the use of R determined 
from polarization transfer experiments, Eq. (|5|) can be expressed as: 

a R = G 2 M {Q 2 )(l+ £ -R 2 ^ +F(Q 2 ,e). (8) 

In other words, in Eq. ([8]) the input parameters are now Gm an d F with R given by Eq. 
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Table 1: The values of the form factors and coefficient A obtained in the fit II, Eq. ([9]), with 
B(Q 2 ) assumed to be 0. 
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Table 2: The values of the form factors and coefficient B obtained in the fit II, Eq. ([9]), with 
A(Q 2 ) assumed to be 0. 

In contrast to Eq. ([7]) which is singular at y = (e = 1), we will assume F(Q 2 ,y) is 
analytic around y = 0. The most natural choice of the parametrization of F(y) will then be a 
polynomial of y as would be obtained in a Taylor expansion. In addition, y n ,n > 1 increases 
when e decreases, vanishes at s = 1 (y = 0) and stays finite when e = (y = 1). This feature 
agrees with the results of the model calculations [9j [TO] . Therefore one may parameterize Eq. 
([8]) as a function of three parameters, Gm, A, and B: 

or = G 2 M (Q 2 ) (l + £ -R 2 ^j + A{Q 2 )y + B(Q 2 )y 3 . (9) 
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Table 3: The values of the form factors obtained with the Rosenbluth fit (fit I). The x 2 values 
are also presented. 
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Table 4: The values of the form factors and coefficients A and B obtained in the fit II, Eq. 
([9]). The x 2 values of this fits are also presented. 

It is natural to expect that both A(Q 2 ) and B(Q 2 ) are smooth functions of Q 2 . Tables [T] 
and [2] show the results of the fits where either A or B is set to be zero. The x 2 ' s obtained 
with these parametrizations are as small as those given in the Rosenbluth fit of Table O The 
values of A or B decreases but the ratio between A (or B) and G 2 e /t increases when Q 2 
increases. It is difficult to judge which one gives a better fit because at Q 2 =2.5 and 4.0 GeV 2 
the x 2 with B = is smaller, while for Q 2 = 1.75, 3.25, and 5.0 GeV 2 the fit with A = 
gives smaller x 2 ■ Therefore it is reasonable to keep both A and B in Eq. ([9]). If we further 
parametrize A(Q 2 ) and B(Q 2 ) by 

A(Q 2 ) = aG 2 D (Q 2 ), B(Q 2 ) = f3G 2 D (Q 2 ), G D = - Q l /Q n)2 , (10) 

then the number of parameters reduces from fifteen to seven (the values of Gm at five different 
Q 2 values and two constants a,P). The result of such a fit is presented in Table H] where 
a = -0.221 and (3 = -0.28. 

It is also possible to parameterize F by functions which are not analytic at y = but 
behave smoothly in the region of < y < 1 (0 < e < 1). The logarithmic or double- 
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Table 5: The values of the form factors and coefficients A and B obtained within the fit III, 
Eq. pip . The x 2 values of this fits are also presented. 
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logarithmic functions which are not analytic at y = 0, are such types of functions often 
appear in the loop diagrams pU]. So yln|y| or y(ln|y|) 2 should also be possible choices for 
the TPE contribution. However, the logarithm type functions will be zero when y approaches 
one. This feature does not agree with the model calculations |10j . Therefore we assume that 
F should be parameterized by some combinations of a polynomial of y and some logarithm 
type functions of y such as: 

a R = G 2 M (Q 2 ) (l + £ -R 2 ^j + A(Q 2 )y + B(Q 2 )y(ln \y\) 2 . (11) 

If we again parametrize A(Q 2 ) and B(Q 2 ) by 

A(Q 2 ) = aG 2 D (Q 2 ), B(Q 2 ) = /3G 2 D (Q 2 ), G D = —^—^ , (12) 

we obtain results as given in Table [5] where a = —0.614 and $ = —0.205. 

Let us now compare the results of three different fits, the Rosenbluth fits of Eq. ([T]), with 
Gm and Ge as parameters (fit I), polynomial fit of Eq. ([9]) with parameters Gm and A and 
B (fit II), and the polynomial plus logarithmic function fit of Eq. (|11|) with Gm, A, and B 
(fit III) as the parameters. Their results are presented in Tables El [U and [5l respectively. In 
fit II, A and B are parametrized according to Eq. (|10|) and in fit III, A and B are given by 
Eq. (I12p . It can be seen that the resultant values of x 2 m fits II and III are as small as or even 
smaller than the ones of the Rosenbluth fit. Hence, our choices of the parametrization of the 
TPE effect F indeed can explain the ep cross section data with the values of R = Ge/Gm 
fixed by the polarization transfer experiment data. It is worth emphasizing that there are 
ten input parameters in the fit I but only seven input parameters in the fit II and III. 

Even though the ansatz used in the fits II and III are quite different, there are some com- 
mon features. First, the effect of the I7027 interference is destructive, namely, contribution 
of F is always negative. Consequently, the values of Gm obtained in fits II and III are always 
larger than those of the Rosenbluth fit. In fits II and III, the values of Gm are enhanced by 
about 2% which is very small but nevertheless comparable with the experimental uncertainty. 
In general the values of Gm m the fit III are larger than the ones in the fit II. In addition, 
the values of Ge at higher Q 2 > 2 GeV 2 are much smaller than the results of the Rosenbluth 
fit because the slope of <jr{e) receives large contribution from the TPE effect, as will be 
explained in more details later. 

Fig. Q] shows (Tr{e) at Q 2 = 2.5, 3.25, 4.0, and 5.0 GeV 2 , respectively, for the data from 
Ref. [7]. The dashed, solid, and dotted lines, correspond to the results of the fits I, II, and 
III, respectively. The dashed lines are straight lines and the dotted lines and the solid lines 
have small curvatures. One sees that both of the solid and dotted lines, in the range of 
0.2 < e < 0.95, are close to the dashed lines and both of them show very little nonlinearity 
which agrees with the previous analyses [15\ [TB] . It is interesting to see that even though the 
I7C5 27 interference term F is in principle not a linear function of e, it contributes significantly 
to the slope but very little to the curvature in fits II and III. Furthermore, we observe that 
outside the current data region, the dotted lines show larger nonlinearity than the solid lines 
at very high e region (e > 0.95) . Hence, if the logarithm type functions play important role 
then the nonlinearity of an will become more pronounced at e > 0.95. On the other hand, 
in the region of e < 0.2, the solid lines show larger nonlinearity than the dotted lines. The 
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nonlinearity of both solid and dotted lines are too small as compared to the uncertainty of 
the current data. 

To estimate the contribution of the TPE effect to the slope of o\r(s), one can expand Eqs. 
flU and (HU) around e = 1/2 as in [IS]: 




.1) +C2 ( E _1) 2 + 0((E _1)3 ) , (13) 

where Co, Ci, and C% arise from the interference term F. In fit II, one has 

A + W 4(A + B) 16B 

C0 = ^-' Cl = _ 3Vl~' ° 2 = ^ (14) 

while in the fit III, one has 

4i + ln 2 3B 4i + (4 In 3 -In 2 3)B (16-81n3)S 

C0 " 4^3 ' Cl " 37! ' 2 " 9^ • (15) 

The value of Cq is related to the contribution of the TPE effect to the intercept of or. C\ 
represents the TPE contribution to the slope of or. C2 reflects the curvature of or. Fig. [5] 
shows the value of C\/ {G 2 e /t) as a function of Q 2 , which represents the ratio between the 
contributions of the TPE and the OPE effects to the slope of or, for both fits II and III. One 
sees that the value of C\/(G 2 e /t) increases from a value of ~ 0.3 at Q 2 = 1.75 GeV 2 to about 
4 — 6 at Q 2 = 5 GeV 2 , a huge effect. On the other hand, the value of Cq/G 2 m remains small. 
We stress that it is Ge term, not Gm term one should look for the TPE effect because the 
contribution of the TPE effect to the cross section is only few percents but it dramatically 
changes the slope of the Rosenbluth plots at high Q 2 . If one neglects the TPE effect, the 
values of Ge will be greatly overestimated at large Q 2 . It is interesting to observe that A (A) 
never appears in C2, the coefficient of (e— 1/2) 2 . In [J5], the following form has been adopted 
in their analysis, 



o-r = Po 



l + Pi(e- l -) + P2{e-\f 



(16) 



and they obtained a very small P%. It is not difficult to understand this result because P2 
corresponds to C2/C0 which is proportional to the ratio between B (B) and G 2 M . In addition, 
one notes that if B (B) = 0, then P2 is identically equal to zero. In that case one needs to 
expand or up to 0(e 3 ). The values of P2 are around a few percent as seen in Tables 0] and 
which agree with the findings in [15^ [To] . In general the values of P2 obtained in fit II are 
about three times as large as those of the fit III and they increase slowly with Q 2 . 

For positron-proton scattering the interference term between the one-photon and the 
two-photon exchange amplitudes changes sign and yields a ratio i? ± = o(e + p)/o(e~p) ~ 1 + 
ARe{A 2 ^ /A lr ), where A 2 ^ and A 1 "' are the two-photon and one-photon exchange amplitudes, 
respectively. Owing to the low luminosity of the secondary positron beams, almost all the 
comparisons were made at very low Q 2 , or at small scattering angels, corresponding to e > 0.7. 
New positron measurements [18] will be performed in the near future to measure the TPE 
effect at larger angels and moderate Q 2 values. Here we present our in Fig. [3 for fits 
II and III, where dotted, dashed, and solid lines correspond to results of Q 2 = 1.75,3.25, 
and 5 GeV 2 , respectively. Indeed at low Q 2 and high e region, the R is smaller than a few 
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Figure 1: The reduced cross section ctr at Q 2 = 1.75 GeV 2 . The dashed line is the Rosenbluth 
fit. The solid line represents the fit II. The dotted line represents the fit III. 
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Figure 2: The values of C\/(-G\) in the fit II and III. The solid line represents the fit II and 
the dashed line represents the fit III. 
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Figure 3: The ratio R of the positron-proton and electron-proton scattering cross section 
according to the fit II (left panel) and fit III (right panel). The dotted, dashed, and solid 
curves correspond to results of Q 2 = 1.75,3.25, and 5.0 GeV 2 , respectively. 
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percent. The values of R are not sensitive to the values of Q 2 , for example, in fit II, the 
value corresponding to Q 2 = 1.75 GeV 2 is about 2% lower than the one at Q 2 = 5 GeV 2 . The 
situation is very different at lower e region. The value of R^ would reach more than 10% at 
moderate Q 2 > 1.75 GeV 2 , for both fits II and III, in the region e < 0.4. Consequently the 
new positron measurements |18j will very likely shed new light of the TPE effect in the ep 
scattering if the experimental uncertainty can be reduced with 5% and the e can reach 0.4 
or below. 

In summary, we re-analyze the data of the elastic electron proton scattering to look 
for model-independent evidence of the two-photon-exchange (TPE) effect. In contrast to 
previous analysis of [15] , we parametrize the TPE effect in forms which are free of kinematical- 
singularity, in addition to being consistent with the constraint imposed by crossing symmetry 
and the charge conjugation. Moreover, we fix the value of R = Ge/Gm as determined from 
the data of the polarization transfer experiment. We find that the TPE effect in the slope 
of cr^(e), starting at high Q 2 ~ 2 GeV 2 values, is large and comparable with that coming 
from Ge in the elastic cross section of ep scattering. It becomes dominant when Q 2 reaches 
> A GeV 2 . In addition, the TPE effect is found to behave quasi- linearly in the region of 
the current data, i.e., 0.2 < e < 0.95. It explains why the previous fits do not find clear 
evidence for the nonlinearity in the data of cjr{e). Hence the fact that the current elastic ep 
cross section data shows little nonlinearity with respect to e can not be used to exclude the 
presence of the TPE effect. 

Our values of Ge is much smaller than the results of the Rosenbluth fit. On the other 
hand, our values of Gm are about 2% larger than the results of the Rosenbluth fit. The 
ratio of the positron-proton and the electron-proton scattering cross section is also estimated 
according to our parametrizations. We find that the value of R^ would reach more than 10% 
at moderate Q 2 > 1.75 GeV 2 in the region e < 0.4 which would be of interest for future 
experimental planning. More precise data taken outside the region of the current data e < 0.2 
or e > 0.95 will be required to distinguish the different parameterizations of the TPE effect. 
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